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ABSTRACT 



Frequency domain methods are employed for the design and 
analysis of anti-roll tank systems. A waveslope spectrum of 
the ocean is developed based upon the trochoidal model of 
ocean wares and the Pierson-Moskowitz ocean wave power 
spectra. The response of a typical ship to the waveslope 
spectrum is studied with the aid of a numerical example. 

FORTRAN subroutines for the design of passive and 
pump-activated O-tube anti-roll tanks are developed, and the 
effects of a tank installation on the ship parameters is 
observed. Many controller design methods are investigated, 
and final design is accomplished by exhaustive search of 
closed-loop frequency responses. Frequency responses of 
passive and activated tank stabilized ships are compared 
with that of an unmodified ship to illustrate the 
effectiveness of design. The effects of vessel displacement 
variations due to consumption of fuel and stores on tank 
effectiveness and controller design are discussed. 



4 



TABLE OF CONTENTS 



I. FORMULATION OF THE PROBLEM 8 

A. WAVESLOPE SPECTRUM OF THE OCEAN 8 

B. SHIP ROLLING IN WAVES 12 

1. Definitions and Concepts of Naval 

Architecture 12 

a. Displacement 12 

b. Center of Gravity 14 

c. Center of Buoyancy 14 

d. Metacenter 15 

e. Metacentric Height 15 

f. Righting Arm 15 

g. Righting Moment 15 

h. Metacentric Height and Stability 15 

2. Equations of Motion 16 

a. Natural Period of Roll 16 

b. Roll by Wave Action 17 

c. Equations of Motion 17 

3. Frequency Response 19 

II. SHI? STABILIZATION BY PASSIVE ANTI-ROLL TANKS 23 

A. DEVELOPMENT OF THE EQUATIONS OF MOTION 23 

1. Double Pendulum Analogy 23 

2. Nomenclature 24 

a. Differential Equations 24 

b. Tank Geometry 25 

3. Euler-LaGrange Equations of Motion 27 

4. Calculation of the Energy Functions 28 

5. The Equations of Motion 29 

6. Definitions of Parameters 30 

B. ANALYSIS OF THE EQUATIONS OF MOTION 31 

1. The Self Terms 31 

2. The Mutual Terms 31 



5 



C. ANTI-ROLL TANK DESIGN CRITERIA 32 

D. FREQUENCY RESPONSE OF PASSIVE SHIP-TANK SYSTEM 35 

III. SHI? STABILIZATION BY ACTIVATED ANTI-ROLL TANKS... 40 

A. MODIFIED EQUATIONS OF MOTION 40 

1. Nomenclature 40 

2. Modification of Parameters 41 

3. Equations of Motion 41 

B. CONTROLLER DESIGN AND FREQUENCY RESPONSE 41 

IV. EFFECTS OF VESSEL DISPLACEMENT VARIATIONS 50 

V. CONCLUSIONS 52 

APPENDIX A : Program Listing for Frequency Response of 

Unmodified Ship# with FORTRAN Subroutine 

PARAM 54 

APPENDIX B : Program Listing for Frequency Response of 
Ship with Passive Tanks, with FORTRAN 

Subroutine CALC 58 

APPENDIX C : Program Listing for Frequency Response of 
Activated Ship/Tank System, with FORTRAN 

Subroutine CALCA 63 

APPENDIX D : Glossary for Frequency Response Programs... 68 

APPENDIX E : Waveslope Spectrum of the Ocean 71 

BIBLIOGRAPHY 81 

INITIAL DISTRIBUTION LIST 82 



6 



LIST OF FIGURES 



Figure Page 

1-1 Definitions and Concepts of Naval Architecture . . 13 

1- 2 Frequency Response of Typical Ship 21 

2- 1 The Dependent Variables and Nomenclature 

of Tank Geometry 26 

2- 2 Frequency Response of Ship with Passive Tanks . .39 

3- 1 Block Diagram Development for Activated 

Ship/Tank System 45 

3-2 Frequency Response of Activated Ship/Tank 
System, Proportional Plus Rate Plus 

Acceleration Control 49 



7 



I. 



FORMULATION OF THE PROBLEM 



A ship on the ocean is in constant motion. Ideally, the 
only motions would be translation in the forward direction 
and rotation about the vertical axis, associated with 
propelling and maneuvering the ship. Realistically, there 
are perturbation motions in addition to the desired motions 
in the form of translations along and rotations about the 
longitudinal, transverse, and vertical axes of the ship. 
The purpose of this thesis was to study the rotation about 
the longitudinal axis, known as roll, and arrive at a scheme 
for controlling it. A procedure for designing the 
installation of a U-tube anti-roll tank system was 
consolidated from several sources ^Chadwick, 1950*1 ^Chadwick 
and Klotter, 1 953^j [^Blagoveshchensky, 1 962j j Vachanaratana , 
197 3^] |_Ziegler , 1975"] and presented in the form of a FORTRAN 

controller design methods were 



subroutine. 



Many 



investigated, and the final design was accomplished by 
exhaustive search of closed-loop frequency responses. 



A. WAVES LOPE SPECTRUM OF THE OCEAN. 

Before a meaningful analysis of ship motions in water 
could be conducted, the motions of the water itself had to 
be understood. Because a ship tends to ride normal to the 
surface of the water, the most significant physical property 
of waves with respect to ship motions is the waveslope, ^ . 
At the crest and in the trough of a wave, the waveslope is 
zero, while at some point nearly midway between those 
extremes, the waveslope is maximum. 

Several theories of wave motion have been advanced, 
including the simple harmonic wave; the Stoke's Wave, which 
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is a harnonic wave with a second harmonic added; and the 
trochoidal wave. Geometrically, a trochoidal curve is 
constructed by rolling a circle of radius R under a 
horizontal line, and tracing the locus of a point a radial 
distance c<R from the center of the circle. The trochoidal 
curve may be steep or shallow depending on the ratio of r:R, 
and corresponds very closely to actual observed ocean waves, 
where 2 ir R = L = wavelength, and 2r = h = trough- to-crest 
wave height. 

Ocean waves are primarily caused by the transfer of 
kinetic eiergy from the wind, which is the result of the 
transfer of thermal energy in weather phenomena. As the 
wind blows over a flat surface of water, the water particles 
at the aic-sea interface are excited by frictional forces, 
and tend to move downwind. The freedom of motion of these 
particles is limited by the fluidic properties of water, and 
thus they tend to accumulate in successive ridges, or 
ripples, on the sea surface. Once these ridges are formed, 
further wave buildup continues by a different mechanism. 
On the upwind side of a ridge, a small local high pressure 
region develops which exerts a downward force on the water 
surface, forming the wave trough. On the downwind side of a 
ridge, a corresponding small local low pressure region 
develops which draws the water surface upward, forming the 
wave crest. The wave continues to grow as long as the 
wind continues to supply energy in excess of that consumed 
by the internal friction of the water in the wave, which it 
does until the wave veLocity equals the wind velocity. 

Three key factors govern the development of wave 
systems. They are the initial difference between the 
velocities of the wind and any existing wave systems; the 
length of time (duration) that the wind continues at a given 
velocity; and the fetch, or distance of relatively uniform 
depth water over which the wave will run before being 
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altered by a change of depth or destroyed by encountering a 
land mass. (The passage of a vessel through an ocean wave 
system does little to alter the growth or shape of the 
system . ) 

The rate of increasing wave height is a function of the 
initial difference between the velocities of the wind and 
the existing wave system. A sudden change of wind velocity, 
such as that associated with the passage of a front, can 
cause the waves to increase at a rate of one or two feet per 
minute. There is a physical limit to the increase in wave 
height, aid it is defined in terms of the wave height to 
wave length ratio, h:L=1:7. Near this ratio limit, excess 
wind enecgy blows off the tops of the waves, forming 
turbulence known as "white caps." 

The rate of increasing wave length generally lags behind 
the rate of increasing wave height. Thus in regions of 
short fetch, a sudden change of wind velocity may never 
produce waves with long wavelengths due to the limiting 
ratio. In contrast, on the open ocean where there is a long 
fetch, a change of wind velocity can produce long waves of 
relatively low height. For waves of length L>100 feet, there 
is a diminishing likelihood of ever approaching the maximum 
value of the ratio. In practice, for wavelengths between 
500 and 500 feet, the ratio rarely exceeds 1:20; for the 
longer ocean waves (of length over 1000 feet) , the ratio 
drops to 1:50. 

The above discussion of wave development only hints at 
the complexity of attempting to define the waveslope 
spectrum of the ocean. Conferences have been held, and many 
books have been written on the subject of ocean spectra. 
_Conference, Easton, Maryland, 196lJ (^Inoue, 1 966^] 
Comstock, 1 967j |^Gillaer, 196 9^j ^Ziegler, 197 5^j Many ocean 
spectra plot relative wave energy as a function of frequency 
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with "sea state" as a parameter. Others define wind 
velocity and duration as parameters. Still others define 
"significant wave height" as a parameter, with many 
definitions of "significant wave height" itself. Spectra 
have been constructed for given locations citing mean wave 
height, or mean squared wave height, or RMS wave height as 
parameters . 

Two other factors which affect the effective ocean 
spectra are vessel speed and angle of encounter. To 
illustrate these effects on the ocean spectrum, consider a 
ship proceeding with velocity identical to that of a 
single-frequency wave train. If the ship heads directly 
into the approaching wave train, the apparent frequency of 
the wave doubles, but the roll coupling with the ship goes 
to zero. When the ship travels parallel to the wave train, 
which is to say the ship is in a beam sea, the apparent 
frequency of the wave is the same as the actual frequency of 
the wave, and the roll coupling is maximum. When the ship 
heads in the same direction as the wave train ("following 
seas" conlition) , the apparent frequency of the wave drops 
to zero, as does the roll coupling. If B denotes the ship's 

O 

angle or encounter with the wave train, where B=0 for a 
following sea, B=90° for a beam sea, and B=180°for a head 
sea, then the apparent waveslope is related to the actual 
waveslope, and the apparent frequency of the sea is related 
to the actual frequency by the expressions j^Ziegler, 1975^: 



'P (B) = ^ sin B 

max 

(B) = w (1 . - w u /g) cos B 

sea sea ship 

u = forward velocity of the ship 

ship 

g = acceleration due to gravity 

Since the ultimate reason for this study was 



(D 

( 2 ) 



to aid in 



11 



the design of a general-purpose anti-roll stabilizer, the 
ocean spectrum chosen was a wave height power density 
spectrum, with sea state and significant wave height (1/3 
highest waves) as parameters. The waveslope spectrum was 
than computed on the basis of the trochoidal model of the 
ocean, where 



1/f = 


2 TT / W 
2 


seconds 


(3) 


5.118 


T 


feet 


(4) 


(h/L) 


x 57.3 


degrees 


(first approximationj 



B. SHIP ROLLING IN WAVES. 

A ship on the ocean moves, in general, with six degrees 
of freedom relative to the coordinate axis system. By 
convention, the positive x-direction is forward along the 
longitudinal centerline of the ship, the positive 
y-directioQ is horizontally outward on the starboard beam, 
and the positive z-direction is vertically downward. The 
center of the coordinate axis system is located at the 
vessel center of gravity. Rotation about the x-axis is 
known as roll, with positive roll angle, <t> , defined as an 
inclination to starboard. Roll can be considered separately 
from the other degrees of freedom, which reduces the study 
of ship motions to a two-dimensional problem with one degree 
of freedom. 

1 . Definitions and Concepts of Naval Architecture . 

The following definitions and concepts of naval 
architecture ^Comstock, 1967^ J^Gillmer, 1969^j are given for 
the benefit of the reader, and are illustrated in Figure 

1 - 1 . 

a. Displacement. 

The displacement of a vessel, A , is defined as 
the weight of the volume of water displaced by the vessel. 
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t 




G = Center of Gravity 
B = Center of Buoyancy 
M = Metacenter 
GM = Metacenter Height 
GZ = Righting Arm 

<f> = Ship Roll Angle Relative to Earth Reference 

WL = Water Line 



Figure 1-1. Definitions and Concepts of Naval Architecture 
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Due to the Principle of Archimedes, the displacement equals 
the weight of the vessel itself. Displacement is always 
given in long tons, where one long ton equals 2240 U.S. 
pounds. The weight of the water displaced equals the weight 
of the vessel, independent of the density of the water. The 
volume of water displaced, however, is a function of the 
density of the water. For purposes of most calculations, 
the density of fresh water is given as 36 cubic feet per 
ton, whereas the density of salt water is given as 35 cubic 
feet per ton. 

b. Center of Gravity. 

The center of gravity, G, is defined as the 
point at which the summation of the moments of all the 
weights in a vessel with reference to any axis passing 
through that point equals zero. The center of gravity of 
most warships of usual form is generally located near the 
designer's waterline, in or near the plane of the midship 
section. For purposes of buoyancy and stability analysis, 
the entire weight of a vessel is considered to be 
concentrated at, and the force of the weight acts vertically 
downward through, the center of gravity. 

c. Center of Buoyancy. 

The center of buoyancy, B, is defined as the 
geometric center of the underwater form of the vessel. The 
force of buoyancy, equal in magnitude to the force of the 
displacement of the vessel, acts vertically upward through 
the center of buoyancy. When a vessel is in static 
equilibrium, the center of buoyancy is located directly 
under the center of gravity, such that the force vectors of 
displacement and buoyancy are equal and opposite. Due to 
this definition of the center of buoyancy, as the condition 
of a vessel varies from that of static equilibrium, the 
location of B varies also. 
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d. Metacenter. 

As a ship inclines, the shape of the underwater 
form changes, and the center of buoyancy migrates to remain 
at the geometric center. The migrations form a locus which 
is dependent upon and directly related to the shape of the 
ship's hull. This locus is generally elliptic and as such 
cannot have a fixed center, but for an infinitesimal 
increment of the locus the instantaneous center, called a 
metacentec, can be located. The locus of metacenters is 
geometrically an evolute, and is called the metacentric. 
For small angles of inclination, up to 10° for most ships, 
the metacenters are co-located at a point known as the 
Metacentec, M. 

e. Metacentric Height. 

When a ship is upright, that is, when the angle 
of inclination <t> = 0° , the centers M, G, and B lie in a 
straight vertical line over the keel. The distance between 
G and M is known as the metacentric height, GM. If M is 
above G, GM is defined as being positive. 

f. Righting Arm. 

When a ship is inclined, the center of buoyancy 
is no longer on the original vertical centerline. A point Z 
is described by the intersection of a line drawn vertically 
upward from B with a line drawn horizontally away from G. 
The distance GZ is known as the righting arm, and for small 
angles of inclination <*> , 

GZ = GM sin <t> = GM 4> (6) 

g. Righting moment. 

The righting moment, U r , is the product of the 
force (displacement) and the righting arm: 

M r = A GZ = A GM 4> (7) 

h. Metacentric Height and Stability. 

The metacentric height is a measure of the 
statical stability of a vessel. A large positive GM 
corresponds with a pair of high-frequency complex conjugate 
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poles in the left half of the s-plane, far removed from the 
jw-axis. A small positive GM corresponds with a pair of 
low-frequency complex conjugate poles in the left half of 
the s-plane, close to the j« -axis. If GM = 0, then there 

are two poles at the origin of the s-plane. A negative GM 

puts poles in the right half of the s-plane. 

2 . Equations of Motion . 

a. Natural Period of Roll. 

In still water, a ship can be made to roll by 
applying and removing external heeling moments. In the 
inclined position, the righting moment is equal and opposite 
to the heeling moment, and the potential energy of the ship 
is equal to the work done by the external force causing the 
inclination. When the external moment is removed, the 

righting moment produces rotation of the ship towards the 
upright position. The potential energy is converted to 

kinetic energy such that when the ship is in the upright 
position, assuming unresisted roll, all the potential energy 
has been converted to kinetic energy. As a result, the ship 
will continue to rotate past the upright position until all 
kinetic energy has been reconverted to potential energy, at 
the opposite limit of roll. Assuming no loss due to 

friction, the ship would continue to oscillate, or roll, 
indefinitely with constant amplitude. The period of these 
oscillations can be found from simple harmonic motion 
relationships to be (^Comstock, 1967^ 

T^ = (2"-k) 4/gGM" 

= (1 .108k) / v/gF 

« (0.44b) /J GM (8) 

where 

2 

g = acceleration due to gravity, ft/sec 
k = radius of gyration of mass of ship, feet 
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b = extreme beam of ship, feet 



b. Roll by Wave Action. 

The cause of ship rolling in a seaway is 

primarily the unbalanced moments which result from the 

center of buoyancy shifting in response to the change of the 
waterplane, and thus the shape of the underwater portion of 
the hull, as waves pass under the ship. A ship will usually 
attempt to ride normal to the surface of a wave. A ship 
also moves in a circular orbit as do the particles of water 
in a wave. The orbits induce centrifugal forces, which are 
opposed by the dynamic forces in water. 

Rolling in waves is a forced oscillation, and the ship 
will tend to oscillate at or near the frequency of the 

waves. If the waves are not uniform, which is usually the 
case, a ship will attempt to oscillate with its natural 

frequency. The overall roll, therefore, is a combination of 
the ship’s natural and the wave's periods, with the wave 
influence being generally more predominant. 

c. Equations of Sotion. 

A ship can be modeled as a second-order, 
spring-mass-damper system. ^Chadwick, 1950^J ^Comstock, 1967J 
The spring constant, or static righting coefficient, was 
discussed earlier in this paper in the definition of the 
righting moment. The static righting coefficient is 

K = A G M (9) 

s 

The inertia term, or roll inertia coefficient, is the 
product of the mass of the ship and the radius of gyration 
of the mass of the ship squared: 

2 

J = ( A/g) k (10) 

s 

The radius of gyration is found by timing the natual 
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period of roll in protected (still) water, and massaging the 
equation for roll period as follows: 



k = (T^^/gSH ) /2 tt 

= (T^T^T )/1. 108 (11) 

The coll damping coefficient is derived from the 
observation that the damping ratio of a typical ship is 
approximately 

r = 0.1/v = 0.0318 (12) 

ship 

3y working backwards through the second-order system 
relationships, the roll damping coefficient was found to be 



2 

B = 0.1 (T. AGS ) /» (13) 

s * 

On the forcing function side of the equation of motion, 
following the classical usage of Froude and Krylov, 
^Chadwick and Klotter, 1953^J the ship-sea coupling 
coefficient was considered to be approximately equal to the 
static righting coefficient. Therefore, the equation of 
motion for a ship in a seaway is: 



J £ + B + K <t> =K'/' (14) 

s s s s 

Dividing through by the roll inertia coefficient and 
performing the LaPlace Transform operation, the system 
transfer function was obtained: 



[• 

[s 



2 

+ 

2 



+ 



(B /J ) S + (K /J ) <*> (3) = (K /J ) *(S) 

S S S S -J s s 

2 2 

(2 r u ) s + w "1 4> (s) = « * (s) 

s s s J s 



(15) 

(16) 



4> (s) 



* (5) 



2 

s 



2 

U! 

3 

* (2 f u ) s ♦ 
s s 



~1 

U) 

s 
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(K /J ) 



(17) 



5 s 



2 

S + (B /J ) S ♦ (K /J ) 



S S S -ft. 



It was useful in the analysis to normalize the equations 
such that the ship*s natural frequency was unity. This was 
accomplished in the following manner: 

Starting with the general equation of motion, divide 
through by the static righting coefficient: 



s s s 

The LaPlace Transform operation was performed, and the 
system transfer function (normalized) was obtained: 



3 . Frequency Response. 

In order to study the frequency response of a vessel 
to the waveslope spectrum of the ocean, it was necessary to 
introduce a numerical example. The example chosen was a 



(J /K ) <t> + (B /K ) <t> + <p = 'P 



( 18 ) 



s s s s 



Let t' = cj t = (/K /J ) t 

s s s 

Def ine d <t> /dt = <t> 



d <t> /dt ' = <t>' 
2 2 

d <t> /dt ' = 4>“ 



Then 4> = (Jk /J ) <t>‘ 



s s 



’<P = (K /J ) <*>“ 



s s 

The equation was then rewritten as 
<t>“ + ( B / K J ) <*>' + <t> = t 



( 19 ) 




♦ (B / K J ) s + <|> (s) = * (s) 



s s s 




( 20 ) 



» (s) 
* (s) 



2 

s ♦ (B / K J ) s + 1 



1 



( 21 ) 



s s s 
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well used example l^3lagoveshchensky , 1962J, which was: 

L = 49.1 meters (161.0 feet) 

b = 9.0 meters ( 29.5 feet) 

A = 936.0 metric tons (921.2 long tons) 

GM = 0.73 meter ( 2.4 feet) 

I. = 3.63 seconds 

9 

The radius of gyration of the mass of the ship, and the 
coefficients of the equation of motion were computed as 
follows : 

k = (8. 637(9. 81) (0.73) ) /2* = 3.68 meters 

2 2 
J = (936/9. 81 ) (3 . 68) = 1289.021 tons-meters-sec 

s 

2 

B = (0. 1) (8.63) (936) (0.73)/* =59.7 tons-meters-sec 

3 

K = (936) (0.73) = 683. 28 tons-meters 
s 

These values were then substituted into the transfer 
function, the roots of the characteristic equation were 
found, and the frequency response was determined. 



* (s) 

'i' (s) 



2 

s 



0.53 

+ 0.046s + 0.53 



( 22 ) 



roots (poles) at s = -0.02 + j 0.73 

natural frequency w = 0.73 

s 

damping ratio C = 0.032 

s 

The frequency response is graphically illustrated in 
Figure 1-2. The Digital Simulation Language (DSL/360) 
program which generated the plot is listed in Appendix A. 
The zero decibel low frequency response demonstrates the 
natural tendency of a ship to ride normal to the surface of 
a wave. The resonant peak which occurs at 0.724 rad/sec has 
a gain of 23.9 db, or a magnification factor of 15.6. 
According to the waveslope spectrum of the ocean, tabulated 
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Figure 1-2. Frequency Response of Typica Ship 



21 



in Appendix E , for a sea of 15 foot significant wave height 
(sea state 6) the maximum waveslope angle at 0.724 rad/sec 
would be 2.23 degrees. Therefore, the ship would experience 
a maximum roll angle of 34.8 degrees, relative to the 
horizon. This severity of roll causes crew discomfort, and 
limits the ability of the ship to conduct operations, such 
as gunfire support, flight operations, or small boat 
operations . 
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II. SHIP STABILIZATION 3Y PASSIVE ANTI-ROLL TANKS 



A. DEVELOPMENT OF THE EQUATIONS OF MOTION 

In ordar to apply servo theory to the problem of ship 
stabilization, the physical system was. assumed to be linear, 
with lumped parameters and constant coefficients. It was 
also assumed that the rolling degree of freedom could be 
considered independently of the other degrees of freedom, 
such that the ship oscillates as if about a fixed center of 
rotation. The problem thus became a two-body problem, the 
ship and the tank-fluid, each with one degree of freedom. 
The tank-fluid was idealized as a one-dimensional filament, 
with provision for varying cross-section but with constant 
velocity across each cross-section. These assumptions 
required small angles of ship roll but allowed for moderate 
angles of tank-fluid level, both being reasonable in a 
stabilized ship. 

1 . Double Pendulum Analogy. 

As discussed earlier in this paper, a ship at sea 
behaves like a damped oscillator, and was modeled as a 
second-order spring-mass-damper system. Another appropriate 
model which was considered was the pendulum model Chadwick, 
195 O^j ^Chadwick and Klotter, 1953""], with a spring constant, 
an inertia constant, and a damping constant. 

The stabilizer chosen was a U-tube anti-rolling tank 
system. The water or other fluid in the U-tube acts as a 
hydraulic pendulum. Since the tank was to be contained 
within and wholly supported by the ship, the ship-tank 
system that resulted was a pendulum supporting a pendulum, 
or a double pendulum. In general, this kind of system leads 
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to two sacond-order differential equations with coupling 
ter as . 

2 . Nomenclature. 

There are two types of nomenclature needed here. 
These are the parameters and variables of the differential 
equations, and the dummy variables and parameters which 
characterize the tank geometry. 

a. Differential Equations. 

The differential equations are second-order in 
nature, and have forcing functions. The variables and 
parameters are as follows: 

Dependent variables of the forcing side: 

= instantaneous effective waveslope 
Parameters of the forcing variables : 

K = ship-sea coupling coefficient 

S3 

Torques produced on the forcing side: 

K ^ = disturbance torque due to waves 

ss 

Dependent variables of the homogeneous side: 

4 > = ship roll angle relative to the horizon 
0 = tank-fluid level angle relative to the ship 
Parameters of the homogeneous variables: 

Ship self-terms ; 

j» = roll inertia coef., unmodified ship (no tanks) 

3 

B’ = damping coef. of unmodified ship 

3 

K' = static righting coef. of unmodified ship 

3 

J = roll inertia coef., modified ship (with tanks) 

3 

3 = damping coef. of modified ship 

3 

K = static righting coef. of modified ship 
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Tank self-terras: 



J = equivalent inertia coef. of tank-fluid 
= skin-friction damping coef. of tank-fluid 
= overall damping coef. of tank-fluid 
= equivalent static coef. of tank-fluid 

Hitual coupling terras: 

J = mutual-coupling inertia coefficient 

3 t 

3 = mutual-coupling damping coefficient 

3 t. 

K = mutual-coupling static coefficient 



It is a feature of the double pendulum that 
K = K 

3 1 t 

B = 0 
. St 

As discussed earlier in this paper 
K = K 

3 S S 



b. 



depends upon 
given below, 
s 



R 



r(s) 



Tank Geometry. 

The notation used in the derivation which 
tank geometry is shown in Figure 2-1, and is 

= a generalized coordinate, the curvilinear 
distance along the fluid trajectory. s=0 at 
the tank-fluid surface in the starboard tank, 
and s=S at the tank-fluid surface in the port 

o 

tank, when 0 =0. 

= lever arm of tank, the horizontal distance 
from the vertical centerline of either side 
tank to the vertical centerline of the ship. 

= the radial distance from the virtual center 
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Figure 2.1. The Dependent Variables and Nomenclature of Tank 
Geometry 
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the fluid 



P 

D(S) 

a ( s) 

7<S) 

* 

a 

a(s) 

g 



of rotation to any point on 
trajectory. 

density of fluid in tanks. 



the pec pendicular 
center of rotation 
trajectory, 
angle between r (s) 
angle between r(s) 
tank cross-section 
tank cross-section 
acceleration of gr 



distance from the 
to tangent at any p 

and D (s) . 
and deck of ship, 
at s=0 and s=S. 
at any point s. 
vity . 



vir 

oin 



tual 
t on 



3 . Suler-LaGrange Equations of Motion. 

Tie two degree of freedom system of ship and 
tank-fluid was best treated by means of the Euler-LaGrange 
equations in generalized coordinates. The generalized 
equation is ^Chadwick and Klotter, 1953^j 

a i a u a p 

+ aq + <3q + a q = Q (23) 

i i i i 




where Q ^ is the so-called g 
in the ith equation, and 
equal to the rate of doing 
non-conservative forces 
displacement. The paramete 
are : 



eneralized non-conservative f 
i=1,2. The generalized fore 
virtual work against the exte 
when g^ is given a vir 
rs of the Euler-LaGrange egua 



orce 
e is 
rnal 
tual 
tion 



T 

0 

D 

f 

0 



2 


mv /2 


= kinetic energy 


(24) 


2 


/f.dx 


= conservative potential 
(due to gravity) 


energy 

(25) 


2 


2 

0 v /2 


= dissipative energy 

(due to skin friction. 


(26) 

etc. ) 



conservative force 
damping parameter 
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4. Calculation of the Energy Functions. 

For each type of energy (kinetic, conservative 
potential, and dissipative) , the contribution of the ship 
alone and the contribution of the tank-fluid alone was 
calculated, then the contributions were added to arrive at 
the total energy. In the regions where the fluid surfaces 
travel, namely the vertical legs of the U-tube, it was 
assumed that the tanks were straight, parallel, and of 
uniform cross-section. This meant that a (s) and D{s) were 
constants in these regions, thus for any expressions which 
depended only on the variables a(s) and D(s), the limits of 
integration were fixed at zero and S, without approximation. 



Calculation of kinetic energy: 
S+R tan d 



fluid 



ship 



total 



f p a fsj ( n? r a ~ j 2 + r <*> 

J L L 

R tan d 
•2 

J' 4> /2 
s 

T + T 

fluid ship 



r (s) 



+ 2 



e r a 

a (s) 



<f> D (s) 



.ds 



(27) 

(28) 
(29) 



Calculation of potential energy: 



a 



fluid 



p g A R J' a^(s) r^|s) sin7 ds 



1 - cos <t> 



ship 



+ p g A R 
2 

K' <t> /2 
s 



(30) 



2 tan ff sin4>- tan Q cos 4> 



(J =0 + U 

total fluid ship 



(31) 

(32) 
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Calculation of dissipative energy: 





.2 

B' 9 /2 


(33) 


fluid 


t 


_ 


•2 

B •<*> /2 


(34) 


ship 


s 






D + D 


(35) 


total 


fluid ship 



5 . The Equations of notion. 

The equations of motion of the ship and tank system 
come from equation (23) , with q^ = 4> , and q^ = 9 . 
Substituting in equations (27) through (35) and simplifying. 



and 


letting 


= K 


j f and Q ? 


= o. 


the equations 


of motion 


are 




Ship : 


(J 


$ + B <> *K 


<*> ) * 


(J 9 


+ K d ) 


= K ^ 


(36) 






s 


s 


s 


St 


t 


s 




Tank: 


(J 


$ *K <t> ) 


+ (J 


+ 

DO 


+ K d ) 


= 0 


(37) 






St t 


t 


t 


t 






Obtaining 


the 


system 


tran 


sf er 


function from 


the 


equations of 


motion was so 


mewhat 


more 


complicated tha 


n it 


was 


for the s 


;hip 


without t 


anks. 


First , 


the 


equat ions 


were 



divided through by the respective inertia coeficients, and 
the LaPlace Transform operation was performed: 

2 2 
Ship: fs + (B /J ) s+ (K /J )1<t>(s) + f(J /J )s + (K /J )le(s) 
u ss ss- 1 u sts ts- 1 

= (K /J ) *(S) (38) 

S s 

2 2 

Tank: f(J /J )s + (K /J )l<t>(s)*rs + (B /J ) s+ (K /J )l©(s) 

“-stt tt- 1 L tt tt- 1 

= 0 (39) 

<t> (s) was isolated by applying Cramer’s Rule, and the 
transfer function was obtained: 
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4 


B J +J B 
St St 


3 


K J +J K +B B -2K J 

St St St t St 


= s + 


r 


s + 


2 




j j -j 




J J -J 




S t St 




_ S t St 



B K +K B 
St St 




r 2 ~ | 

K K -K 
St t_ 


IT 

J J -J 

_ S t St _ 


s + 


2 

J J -J 

_ S t St _ 


2 


/J ) S + (B 


/J 


) s+ (K /J ) 


S S 


t t 


t t 



* (3> A 



( 40 ) 



(41) 



6. Definitions of Parameters. 



The parameters of the equations of motion are given 
in terms of previously defined parameters and dummy 
variables as follows: 



Ship self-terms: 



J = 
s 


S 

J' + par y*[a(s)/A][r(s)/R] ds 
0 


• (42) 


B = 


B* 


(43) 


s 


s 




K = 
s 


s 

K' + pgAR y|a (s) /aJ jj: (s) /R^j sin 7 ds 
0 


(44) 


Tank 


self-terms : 

e 




J = 
t 


V 

PAR J A/a (s) ds 

o 


(45) 


B = 


V 

B* 


(46) 


t 


t 


K = 


2 

2 P g A R 


(47) 
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Mutual-coupling terms: 



J = 
St 


s 

p A R J D (s) /H ds 
0 


(48) 


B = 
3 t 


0 


(49) 


K = 
3 t 


K 

t 


(50) 



B. ANALYSIS OF THE EQUATIONS OF MOTION. 

Having arrived at a set of equations of motion, it was 
necessary to conduct a brief analysis of the equations to 
give them physical significance. 

1 . The Self Terms . 

Looking first at the self-terms of the homogeneous 
equations, it was observed that the ship acts as a simple 
pendulum if the tank-fluid does not move, and that the 
tank-fluid acts as a simple pendulum if the ship is 
prevented from moving. Each of these pendulums has its own 
natural frequency. These two frequencies are important 
factors in the stabilization problem. 

2 . The Mutual Terms. 

There is a third critical frequency associated with 
the ship-tank system, called the secondary resonance 
frequency, which is due to the mutual-coupling terms. The 
mathematical significance of the effect became apparent when 
the ship equation, equation (36) , was rewritten with the 
mutual teems on the right side as follows: 

J £ +B +K <j> = K + - (J 9 +K 9 ) (51) 

s s s s st t 
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This equation shows the ship as a simple pendulum, acted 
upon by a net torque, which is the torque due to waves minus 
the torque due to tank-fluid motion. If the tank-fluid were 
made to oscillate at the secondary resonance frequency, the 
tank-fluid torque would be equal to the wave torque, and 
thus the ship would be influenced by zero net torque. 

C. ANTI-ROLL TANK DESIGN CRITERIA. 

The easiest tank geometry to design and build is a right 
rectangular O-tube, with the vertical legs having constant 
cross-sectional area, and the horizontal cross-leg having 
constant but smaller cross-sectional area. Given this tank 
geometry, 

a(vertical leg) = A (52) 

a (horizontal leg) =. a 

= (A R .<*>')/ (g - h a) ) (53) 

s s 

D(vertical leg) = R (54) 

D(horizontal leg) = D (55) 

The length of the trajectory of the fluid in the tank, 
S, is dependent upon the lever arm of the tank, R, and the 
height of the fluid column in a vertical leg of the U-tube 
when 9 = 0 0 , h, 

S = 2 (R + h) (56) 

For maximum effect, the tank system must span the 
extreme beam of the ship. Therefore the lever arm of the 
tank is half of the beam, b, diminished by the width of a 
vertical leg, w, 

R = (b-w)/2 (57) 

The relationships derived from the equations of motion 
which define the tank dimensions and parameters are: 
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K /J 
t s 



= X 



(58) 






= it) 



2 a 



J' A/a(s) ds 
0 



(59) 



J /J = X 
st t 



= /(D a)/(R A) ds 

0 

= 2 h + 2(D a/A) 



(60) 



The first ratio, x, reflects the free surface effect of 
the tank-fluid, and should be a small number. An acceptable 
value was found to be 0.05 (jfachanar atana , 1973^. 

To be effective as a stabilizer, the anti-roll tank 
system must be "tuned" to the natural frequency of the ship 
in which it is installed. Therefore, 



w 



= CJ 



K /J = K /J 
t t s s 



(61) 



The third ratio, X, is dependent upon the placement of 
the tanks in the ship relative to the center of rotation. 
This ratio can be forced to zero (J^ forced to zero) by 
placing the horizontal leg of the U-tube above the center of 
rotation by the amount 



-0 = h A/a 



(62) 



There are no hard and fast rules for determining a value 
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of h, except that the tanks must fit within the hull of the 
ship, and the fluid in the tanks should not be permitted to 
either contact the top or completely drain out of the bottom 
of a vertical leg of the tank at the maximum expected 
tank-fluid level angle, ^ m3l . 

Since J depends upon tank geometry, the approximation 
must be made that J = J' , for purposes of initiating 
design. Then 

K = X J» 
t s 

2 

= 2 p g A 3 

2 

= P g 1 w (b-w) /2 (63) 

A reasonable value of tank width should be chosen. Then 

2 

1 = 2K^/pgw (b-w) (6U) 

Making the approximation that the natural frequency of 
the modified ship is the same as the natural frequency of 
the ship without tanks, and picking a value for h, equation 
(53) is solved for a, and then equation (62) is tested for 
feasibility. If the value of D so computed cannot be 
physically realized, then a realizable value of D must be 
chosen. This merely means that J S | will not be identically 
zero, and the overall analysis will take more care. 

The tank damping parameter is derived by working 
backwards through the second-order relationships for the 
hydraulic pendulum, noting that the damping ratio for the 
tank-fluid is approximately equal to the damping ratio of 
the (unmodified) ship: 



f - 0 . 1/ x 

B =0.2 JK J ~ / x 
t v t t 



(65) 

( 66 ) 



Once the tank system parameters have been determined, 
the modified ship parameters are found. Over the range s=0 
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9 



l 2 

to s=h, r = R + 
range s=h to s=h+R, 



(D + h-s) 

J 



2 



and 



= (R+h-s) 



2 



r sin 7 = D+h-s; over the 
7 

+ D. , and r sin 7 = D . 



2 

J = J* + 2 P A R (I+I) (67) 

s s 12 



I 

1 



[< 



2 2 2 
3R +2D ) h + 2Dh + 




( 68 ) 



I 

2 


= a (R 


+ 3 D )/( 3AR ) 


(69) 


K 

3 


+ 

- V) 

II 


2 p g A R (I + 1 ) 
3 4 


(70) 


i 

3 


= - (2D 


2 

+ h ) /( 2R ) 


(71) 


T 

a. 


= Da/A 




(72) 



Ihe design steps should then be retraced, removing the 
approximations. Normally, after two or three iterations of 
the design procedure, the calculated values converge on the 
actual values, and a final design has been achieved. The 
FORTRAN subroutine used in this iterative design procedure 
is included in Appendix B to this thesis. 



D. FREQUENCY RESPONSE OF PASSIVE SHIP-TANK SYSTEM. 

The same numerical example used in the frequency 
response study of a vessel on the ocean was used for the 
design example and frequency response study of a ship with 
passive anti-roll tanks to the waveslope spectrum of the 
ocean : 

2 

J 1 = 1289.021 metric tons-met ers-sec 

3 

3* * 59.746 metric tons-meters-sec 

3 

K' = 683.280 metric tons-meters 

3 

The following arbitrary choices were made: 
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w 

h 

X 

P 3 

Then 

K 

t 

1 



A 



1 . 0 meter 
1.3 meter 



0.05 

1.025 me trie 



3 

tons/meter 



(salt water used) 



(0.05) (1289.021) 

64.451 metric tons-meter 

2 

(2) (64. 451) / (1.025) (1) (9-1) 

1.965 meter 



w 1 
1.965 



2 

meter 



a 




2 

= 0.457 mater 



Testing equation (52) for feasibility, 

-D = (1.3) (1 .965) / 0.457 

= 5.59 mater 

Since this tank placement would probably not be possible 
in the ship being studied, a more realistic value, still 
above the center of rotation, was arbitrarily chosen: 



-D = 1.0 meter 



Then X is non-zero, as indicated. 

J /J = X 
3 1 t 

= 2 []l. 3 - (0.457) (1 .965)]] 

= 2.135 



To summarize, the tank dimensions found thus far were 
w = 1.0 mater 

1 = 1.965 meter 
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k 


= 1.965 meter 






2 


a 


= 0.457 meter 


a 


= 4.0 meter 


D 


= -1.0 meter 


Siven 


these tank 


dimensions, the modified 


par araeters 


became: 






= 1289.021 


+ 11.155 


3 




2 




= 1300.176 


metric tons-meters-sec 


K 

3 


= 683.280 


- 3.121 




= 680.159 


metric tons-meter 


k ft er 


four iterations of the design procedure. 


following 


final values 


were obtained: 


w 


= 1.0 


meter 


1 


= 1.982 


meter 






2 


k 


= 1.982 


meter 






2 


a 


= 0.454 


meter 


R 


= 4.0 


meter 


D 


= -1.0 


meter 






2 


J 

3 


= 1300.236 


metric tons-meters-sec 


B 

3 


= 59.746 


metric tons-meters-sec 


K 


= 680.185 


metric tons-meters 


3 


= 124.276 


2 


j 


metric tons-meters-sec 


t 






B 


= 5.722 


metric tons-meters-sec 


t 






K 


= 65.012 


metric tons-meters 


t 




2 


J 


= 266.152 


metric tons-meters-sec 



ship 



3 t 



These values were then substituted into the transfer 
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function, the roots were located, and the frequency response 
was obtained. 

2 

4> (s) 0.523 [ s t 0 ♦ 0460 s + 0.523 ) 

_ -3 3 2 

(s) s +• 0.164 s + 1.485 s + 0.0857 s + 0.441 

(73) 

zeroes at s = -0.0230 + j 0.723 

poles at s = -0.0135 t j 0.642 

and s = -0.0684 ± j 1.031 

The frequency response is graphically illustrated in 
Figure 2-2. The DSL/360 program which generated the plot is 
listed in Appendix 3. The negative decibel low frequency 
response demonstrates the ability of the passive tank system 
to overcome the natural tendency of the ship to ride normal 
to the waveslope, making it ride upright with respect to 
earth reference instead. 

Comparing this response with the response of the 
unmodified ship, given in Figure 1-2, it was observed that 
the effect of adding passive anti-roll tanks was to place a 
pair of complex zeroes at the ship's own natural frequency 
to reduce the resonant peak, from 24 db at 0.7 rad/sec, to 
two lower peaks of 14 db at 0.6 rad/sec and 10 db at 1.0 
rad/sec, and a notch of -11 db at 0.7 rad/sec. In sea state 
6, this translates to maximum ship roll angles of 8.5 
degrees at 0.6 rad/sec, and 15.3 degrees at 1.0 rad/sec, a 
considerable improvement over the 34.8 degree maximum roll 
angle of the unmodified ship. 
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Figure 2.2. Frequency Response of Ship with Passive Tanks 
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III. SHIP STABILIZATION BY ACTIVATED ANTI-ROLL TANKS 



A. MODIFIED EQUATIONS OF MOTION. 

Activating an aati-roll tank system, mathematically, 
only involves adding a forcing function term to the tank 
equation of motion. Physically, however, it can be 
accomplished in basically two different ways. These methods 
are pneumatic and hydraulic. A pneumatic effector operates 
by alternately creating an overpressure or a partial vacuum 
condition in the airspace above the tank-fluid surface in 
the vertical legs of the U-tube, while an hydraulic effector 
generally consists of an arrangement of pumps in the 
horizontal cross duct. The effector type chosen for this 
analysis was a constant speed controllable pitch propellor 
pump, so chosen for its linear blade pitch angle vs. flow 
characteristic [^Blagoveshchensky, 1962^. 

1 . Nomenclature . 

In addition to the nomenclature of the passive tank 
system, the following are the parameters and variables of an 
activated tank system: 

Dependent variables of the forcing side: 

a = instantaneous pump-blade pitch angle (commanded) 

Parameters of the forcing variable : 

B = velocity head-loss coefficient of the pump 
P 

K = static head coefficient of the pump 
o 

Torques produced on the forcing side: 

K a - B 9 = Counter-rolling torque due to pump 

P P 
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2. Modification of Parameters . 

As indicated by the expression for counter- rolling 
torque, the addition of a pump to the tank system has the 
effect of increasing the overall damping coefficient of the 
tank-fluid. All other parameters of the homogeneous side 
remain unchanged. The pump should be chosen such that the 
composite tank damping factor is 0.45 J^Vachanaratana, 1973^: 

(74) 



B = B* + B 

t t p 

B /J = 0.45 
t t 

B = 0.45 J - 3* 

? t t 



(75) 

(76) 



3. Equations of Motion. 

The equations of motion for the activated anti-roll 

tank system come from equation (23) , with Q = K.’/' , Q = K.at 

i— I — 1 Z r 

, q = <t> , and q 9 = 0 [Chadwick and Klotter, 195 3 I . 

Ship: (J d>* + B 4> *■ K <t> ) + (J 0* + K 6 ) = K + (77) 

s s s St t s 



Tank: (J * + Kd>) + (J0+B0 + K0)=Ka (78) 

st t t t t p 

Following the procedure used for the passive ship-tank 



<t> (s) 


was found 


to be 






4> (s) 


= * (s) - 


4> (s) 
2 




(79) 


(s) 


= (K /J ) 
s s 


|— 2 

s + (B /J ) s + (K /J 
L t t t t 


) J ^ (S) 


(80) 






A 






<t> 2 ( s ) 


= (K /J ) 
P t 


f(J /J )s 2 + (K /J f| 
st s t s J 


A(s) 


(81) 


A 


= as defined by ^equation (40) 







B. CONTROLLER DESIGN AND FREQUENCY RESPONSE. 

The pcoblem of ship roll control is a frequency domain 
problem: controlling against a (time domain) stochastic 



41 



disturbance with a generally known (frequency domain) 
spectrum. Therefore the classical approach to control was 
pursued, designing on the basis of feedback compensation. 
It was assumed that roll angle, rate, and acceleration 
signals would be available from shipboard sensors (gyros and 
accelerometers) . The controller transfer function was 
considered to be of the form 



C(s) = k (s) / <f> (s) 



= K s ♦ K s ♦ K 
2 1 0 



(82) 



For purposes of block diagram development and analysis, 
the following block transfer functions were defined, based 
upon equations (79) , (80) , and (81) : 

SS (s) = <j>^ (s) / 'I' (s) 

2 



= K /J fs + (B /J )s + (K /J )1/A 

SS 1 -- tt t t -* 

2 r 2 

s + 2 fw i 
n L - t n 



= C J 



2-j 

s + o> /A 

n -J 



ST (s) = 4> (s) / & (s) 

2 



= K /J |~(J /J ) s + (K /J )!/ a 
pt'-sts ts J 

r 2 2- l 

= K /J Xs ♦ « / A 

p s n- 1 



(83) 



(84) 



SS T (s) = ST (s) / SS (s) 

2 2 2 2 
= K /K F Xs + <J ~]/[s + 2U S + UI 1 (85) 

p s L. n-l L- t n nJ 

G ( 3 ) = K /J fs 2 + (B /J ) s + (K /J )1 

S S L t t t t -I 

A - K /J fs + (B /J ) s + (K /J )1 

SS>- tt tt-1 



w r s 


♦ 


2 r a> s + (Jj 1 


a L 




t n n J 


~ 7~ 


r ^ 


. 2-i 


A — w 

TO 


L s 


+ 2f u s + 

f TO TO — ' 



( 86 ) 
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H(s) 



= SST (s) C ( s) 

4 3 



Where 



= 1/K fxc s +XC s +(XC +wC)s + wCs+cucl 
s'- 2 1 0 n 2 n 1 n 0-1 

2 2 

s + 2 S w s + ^ (87) 

t n n 



= K K 

P 2 

= K K 

P 1 

= K K 

P o 



The "open loop" transfer function, then, was 
F(s) = G (S) / [l + G (s) H ( s)] 



(K /J )fs + (3 /J ) s + (X /J )1 

s s L £ t t t - I 

A (open) 



( 88 ) 



A (open) = A - (K /J )fs 2 + (B /J ) s ♦ (K /J )1 + 

ss 1 - tt tt J 

r- 4 3 2 

+ (1/J ) XC S +X C s + (XC + (K /J ) C ) s 

s L 2 1 0 t t 2 

+ (K /J )C s + ( K / J )C “1 

t t 1 t t 0J (89) 

The overall, or "closed loop" response became 

T(s) = <*> (s) / * (s) 

= F (s) / [l + F(s)] 

r 2 

= (K /J )fs + (B /J )s + (K /J ) | (90) 

s s'- t t t t J 

A (closed) 

I - 4 3 2 

A (closed) * A+ (1/J )|_ XC 2 S + XC ^s + (XC Q + (K^/J^) C^) s 

+ (K /J ) C s + (K /J ) C 1 (91) 

t t 1 t t 0-1 
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The Characteristic Equation 
A (closed) , was 

4 3 2 


of the 


closed loop 


s + A s + A s 
3 2 


♦ As + A 
1 


= 0 

0 








- 


- 


A = 


J J 




B J + J 


B + J C 


3 


s t 

J J +J C 




s t s t st 1 

2 

J J +J J J 




st st 2 




s t st s t 


A 


J J 




K J * J 


K + B B - 2K 


2 


s t sts 

J J +J c 
s t st 2 

J C + K C 
st 0 t 2 

+ 

J J 
s t 


; t s t 

2 

J J +J 
s t st 


A = 


J J 




“B K * K 


B + K C 


1 


s t 

J J +J C 
st st 2 




s t s t t 1 

2 

J J -J J J 

s t st s t 

T 2 - 


A 


J J 




K K - K 


+ K C 


3 


s t 

J J J C 
s t st 2 




s t t tO 

2 

J J -J J J 

s t st s t 



t St 



( 92 ) 



(93) 



(94) 



(95) 



(96) 



The block diagram development is illustrated in Figure 
3-1. The first design approach attempted was the modern 
control theory state variable feedback method. The 
difficulty with this method was the inability to define a 
time-domain objective or cost functional. 

The next design approach attempted was with the Computer 
Aided Design of Optimal Compensators (CADOC) program 
[ Anastassakis , 1977 ]. This method consisted of selecting a 
frequency-domain cost functional and a desired frequency 
response curve, and used the function minimization procedure 
of M.J. Box. There were two difficulties encountered in 
applying this program to the problem of ship anti-roll 
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(a) Direct implementation of equations (79) and (82). 




(b) First state of block diagram reduction. 



Figure 3.1. Block Diagram Development for Activated 
Ship/Tank System 
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(d) Equivalent block diagram for servo theory analysis. 



Figure 3.1. Block Diagram Development for Activated 
Ship/Tank System 
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control. The first problem was that a portion of the 
feedback path in the physical system was fixed by 
architecture/ while the CADOC program required all 
coefficients in the feedback path to be adjustable. The 
second problem involved selecting an attainable frequency 
response for CADOC input. It was concluded that the CADOC 
program was best suited for optimizing a close design, 
rather than for initiating a design. 

Another design approach attempted was in the root locus 
parameter space. Two separate procedures exist here, and 
both were attempted. The first procedure was the standard 
root locus approach, fixing two of the three adjustable 
gains and varying the third. This method amounted to an 
exhaustive search technique, attempting to locate the best 
combination of the three gains to define system roots with 
adequate lamping. 

The second root locus parameter space procedure involved 
using a parameter plane computer program, PARAM A, which 
exists at the Naval Postgraduate School and is available 
from Dr. 3.J. Thaler, but has not been documented in any 
literature to his knowledge. This program varies two gains 
which appear in linear combinations in the coefficients of 
the system characteristic equation, for a fixed value of the 
third gain. The output of this program is a two-dimensional 
parameter plane showing families of curves of constant 
damping ratio and constant natural frequency as functions of 
the two gains. 

The problems associated with both root locus parameter 
space procedures were the exhaustive search nature of a 
three-variable system, and the exclusions of the system 
closed-loop zeroes from the analysis. This was noted to be 
a problem when the gains found by these procedures were used 
in a frequency response study of the activated ship/tank 
system, and large resonant peaks were observed. 
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After attempting all of the above methods and abandoning 
them for the reasons stated, an exhaustive search of 
closed-loop frequency response curves was conducted, using 
the program listed in Appendix C. An acceptable design was 
accomplished when C = 110, C = -150, and C = -80. This 
corresponded to placing the compensator roots at s = -2.4388 
and s = *0.5638, which improved the overall system damping 
and provided the frequency response illustrated in Figure 
3-2. This response was characterized by two low resonant 
peaks, one at 0.582 rad/sec of 2.85 db (linear gain 1.39), 
and one at 1.38 rad/sec of 1.96 db (linear gain 1.25). The 
minimum response between the two resonances occured at 0.861 
rad/sec, and was -4.75 db (linear gain 0.58). 

In sea state 6, the maximum waveslopes at the resonant 
frequencies are 1.44 degrees at 0.582 rad/sec, and 8.11 
degrees at 1.38 rad/sec, and the maximum waveslope at 0.861 
rad/sec is 3.15 degrees. Using this controller design, the 
maximum roll angles that the ship should experience in sea 
state 6 would be 2.0 degrees at 0.582 rad/sec, 1.83 degrees 
at 0.861 rad/sec, and 10.1 degrees at 1.38 rad/sec. The 
wave height power spectral density at those frequencies, 
according to the spectra of Pierson and iloskowitz, are 205 
ft 2 /Hz at 0.582 rad/sec (0.093 Hz), 80 ft 2 /Hz at 0.861 
rad/sec (0.137 Hz), and 12 ft 2 /Hz at 1.38 rad/sec (0.219 
Hz) . Therefore it was concluded that even though the pure 
numbers indicated a poor roll control at the higher 
frequency resonant peak, the waves at that frequency lack 
sufficient power to cause the ship to attain the maximum 
theoretical roll angle. 
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FREQUENCY RESPONSE, ACTIVATED TANKS 
CLOSED LOOP RESPONSE, MAG (DB) VS. LOG W 




Figure 3-2. Frequency Response of Activated Ship/Tank 

System, Proportional Plus Rate Plus Acceleration 
Control 
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IV. EFFECTS OF VESSEL DISPLACEMENT VARIATIONS 



The analysis of ship motions and roll control thus far 
had been conducted under the assumption that the 
displacement of the vessel was constant. Over the short 
term, the assumption may be true enough, but for a ship 
underway, consuming fuel and stores, the assumption cannot 
hold true beyond a few days at best. 

As a ship's displacement decreases, the underwater 
surface area of the hull decreases also, in keeping with the 
Principle of Archimedes. This in turn leads to a lowering 
of the Center of Buoyancy, which results in a lowering of 
the Metacenter. The effect on the metacentric height, 
however, is not so clear, as the exact location of the 
Center of Gravity in the ship depends more on the placement 
of the onboard weight than on the total displacement. In 
general, fuel tanks are located low in a ship; thus 
consumption of fuel is removal of low weight, which shifts 
the Center of Gravity upward, decreasing the metacentric 
height. It should be noted here that ballasting empty fuel 
tanks with sea water has the effect of restoring the ship as 
much as possible to the original condition of statical 
stability which existed at the beginning of the voyage, and 
extends the time frame of applicability of the assumption of 
constant displacement. 

Varying the ship's displacement therefore varies the 
natural frequency of the ship, and effectively de-tunes the 
anti-roll tank system. The only physical adjustment that 
can be made to a passive tank installation underway is to 
vary the quantity of water in the U-tube. This has the 
effect of changing the limits of integration of the design 
equations, and therefore changes the parameters of the tank. 
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Additional adjustments are available on the activated 
tank system. These are the gain of the pump, which is a 
function of the pump rotational speed, and the feedback 
gains from the sensors to the pump controls. Thus, if a 
method exists for determining on-line the various ship 
parameters as they change, the controller for the activated 
tank system can be made adaptive. Various schemes may be 
considered, from Kalman filter type parameter estimators to 
FFT devices, and this field is left as a topic for further 
possible cesearch. 
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V. 



CONCLUSIONS 



In studying ship roll control, it was necessary to 
become familiar with the nature of ocean waves. The 
waveslope spectrum of the ocean which was generated, when 
used in conjunction with the Pierson- Moskowitz wave power 
spectra, gave much insight into the problem. It was learned 
that the sea surface does not and can not assume the slopes 
that an observer on a ship often estimates. 

A typical ship on the ocean magnifies the waveslopes 
near the resonant frequency by a factor of about 15. By 
installing passive anti-roll tanks, this waveslope 
magnification factor was split into two resonant peaks, one 
with a factor of 5, and the other with a factor of 3. 
Activating the tank system further decreased the waveslope 
magnification factors to 1.4 and 1.3, respectively. 

Several assumptions have been made over the years in the 
study of roll control by U-tube tank systems, and some of 
them were sources of problems in the design of physical 
systems. The first assumption was that a ship has constant 
parameters: displacement, metacentric height, and natural 
period. This assumption holds true only if the ship does not 
coasume any fuel or stores, or if these weights are replaced 
with appropriate quantities of suitably placed ballast. 

Another assumption was that the tanks could be placed in 
the ship in such fashion as to drive the mutual -coupling 
inertia coefficient to zero. This assumption had two 
limitations. The first pitfall was the direct connection 
between this assumption and the assumption of constant ship 
parameters. The second pitfall was that the tanks had to be 
placed in an impossible location on the ship to accomplish 
the objective of the assumption. 
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Roll control by pump activated tanks is perhaps the best 
form of roll control for slow moving vessels, or for vessels 
conducting operations either at slow speed or stopped, where 
activated fin stabilizers lose their effectiveness. It 
should be noted that a pump activated tank system can be 
used to cause rolling, as an aid for breaking out of ice. 
It can also be used to correct list, including damage-caused 
list, as a safer alternative to counter- flooding. Therefore 
the pump activated U-tube tank system is the most versatile 
form of coll control that a multi-purpose vessel can have 
installed. 
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APPENDIX A 



Program Listing for Freguency Rasponse of 
Unmodified Ship, with 
FORTRAN subroutine PARAM 



The Fcaguency Response program was written in Digital 



Sim ulation 


Language 


(DSL/360) to 


be 


ran on 


the IBM 


Systera/360 


machine 


at the W. R. 


Church Computer 


Cent er , 


Naval Postgraduate 


School. To use 


the 


program, a 


deck of 


plot identification 


cards must be 


appe 


snded folio 


wing the 


FORTRAN Subroutine, 


in accordance 


with 


the DSL/3 


60 user's 


manual. 












The labels and 


variable names 


used 


in the program are 


defined in 


the Glossary for Fregu 


ency 


Response 


Programs, 


Appendix D 


of this thesis. 
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APPENDIX B 



Program Listing for Frequency Response of 
Ship with Passive Tanks, with 
FORTRAN Subroutine CALC 



The Fre 


guency 


Response program 


was 


written 


in 


Digital 


Simulation 


Language (DSL/360) to 


be 


run on 


the IBM 


System/360 


machine 


at the W. R. 


Church Computer 


Center, 


Naval Postgraduate 


School. To use 


the 


program. 


a 


deck of 


plot identification cards must be 


appended follow! 


ng the 


FORTRAN Subroutine, 


in accordance 


wit h 


the DSL/360 


user' s 


manual. 














The labels and 


variable names 


used 


in the 


program are 


defined in 


the Glossary for Frequency 


Response 


Programs, 


Appendix D 


of this 


thes is. 
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APPENDIX C 



Program Listing for Frequency Response of 
Activated Ship/Tank System, with 
FORTRAN Subroutine CALC A 



The Frequency Response program 
Simulation Language (DSL/360) t 
System/360 machine at the W. R. 
Naval Postgraduate School. To use 
plot identification cards must be 
FORTRAN Subroutine, in accordance 
manual. 

The labels and variable names 
defined in the Glossary for Freq 
Appendix D of this thesis. 



was written in Digital 
o be run on the IBM 
Church Computer Center, 
the program, a deck: of 
appended following the 
with the DSL/360 user's 

used in the program are 
uency Response Programs, 
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TURN 



APPENDIX D 



Glossary for Frequency Response Programs 

MAIN 



SYMBOL 

AO , A 1 ,A2, A3 

BS 

BT 

CLOOP 

COMP 

D 

G 

GIMAG 

GREAL 

H 

IM 

IMH 

JS 

JSI 

JT 

KS 

KT 

LOG W 

MAG 

MAGH 

MDB 

MDBH 

N 

OLOOP 

PHI 



MEANING 

Coefficients of characteristic equation 

Damping coefficient of ship 

Overall damping coefficient of tank-fluid 

Closed-loop transfer function 

Compensator (controller) transfer function 

Denominator polynomial of open-loop transfer 

function 

Open-loop transfer function 

Imaginary part of open-loop transfer function 
Real part of open-loop transfer function 
Closed-loop transfer function 
Imaginary part of open-loop transfer function 
Imaginary part of open-loop transfer function 
Roll inertia coefficient of ship 
Mutual-coupling inertia coefficient 
Equivalent inertia coefficient of tank-fluid 
Static righting coefficient of ship 
Equivalent static coefficient of tank-fluid 
Natural logarithm of radian frequency 
Magnitude of open-loop transfer function 
Magnitude of closed-loop transfer function 
Decibel gain of open-loop transfer function 
Decibel gain of closed-loop transfer function 
Numerator polynomial of transfer function 
Open-loop transfer function 
Phase of open-loop transfer function 



68 



Q 


Denominator polynomial of open-loop transfer 
function 


s 


LaPlace transform variable 


SST 


Ship-tank response in feedback path 





SUBROUTINES 


SYMBOL 


MEANING 


A 


Cross-sectional area of vertical leg of U-tube 


AS 


Cross-sectional area of horizontal leg of 
U-tube 


BEAM 


Transverse dimension of ship (beam) 


BP 


Velocity head-loss coefficient of pump 


BS 


Damping coefficient of ship 


BT 


Overall damping coefficient of tank-fluid 


BT1 


Skin-friction damping coefficient of tank 
fluid 


D 


Perpendicular distance from center of rotation 
of ship to horizontal cross-leg of U-tube 


DABS 


Absolute value (magnitude) of DCAL 


DCAL 


Calculated value of D to force the 


DISP 


mutual-coupling inertia coefficient to zero 
Displacement of ship 


El, E2 


Incremental diffenence between current and 
previously calculated values of ship 
parameters 


GM 


Metacentric height 


GRAV 


Acceleration of gravity 


H 


Height of fluid column in vertical leg of 
U-tube 


11, 12,13,14 


Evaluation of integrals as described in main 
body of thesis 


JS 


Roll inertia coefficient of ship 
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JS 1 
JST 
JT 

KRAD 
KS 
KS 1 
KT 
L 

LEN 
LMD A 
N 
R 

RHO 

RHDG 

TJ 

TK 

TPHI 

WIDE 

WN2 

X 



Roll inertia coefficient of unmodified ship 
Mutual-coupling inertia coefficient 
Equivalent inertia coefficient of tank-fluid 
Radius of gyration of the mass of the ship 
Static righting coefficient of ship 
Static righting coefficient of unmodified ship 
Equivalent static coefficient of tank-fluid 
Longitudinal dimension of a vertical leg of 
the O-tube 

Longitudinal dimension of ship (length) 

Design ratio X 

Humber of iterations of design procedure 

Lever arm of U-tube 

Mass density of tank-fluid 

Weight density of tank-fluid 

Temporary JS 

Temporary KS 

Natural period of roll of ship 

Transverse dimension of vertical leg of O-tube 
Natural frequency of ship-tank system, squared 
Design ratio X 
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APPENDIX E 



Waveslope Spectrum of 



The waveslope spectrum 
the relationships of the 
intended to be used 
Pierson -Moskovitz wave power 
This appendix includes 
program which generated the 
spectra foe sea states 3, 4, 



the Ocean 

of the ocean was generated from 
trochoidal wave model. It is 
in conjunction with the 
spectra . 

the listing of the DSL/360 
spectrum, and listings of the 
5, and 6. 
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RESET LABEL 

LABEL WAVESLCPE SPECTRUM FOR H=15.0 FEET 



WAVESLOPE SPECTRUM OF THE OCEAN 
SIGNIFICANT HEIGHT =3.3 FEET 
******** (SEA STATE 3) ******** 



W 

1. 0000E-01 
1 . 0593E-0 1 

1. 1220E-01 

1 . 1885E-0 1 
1. 2589E-01 
1 .3335E-01 
1. 4125E-01 
1 . 4962E-0 1 
1. 5849E-0 1 
1 .6788E-01 
1. 7783E-01 
1 . 8836E-0 1 
1. 9953E-01 
2. 1135E-01 
2. 2387E-01 
2.37 14E-01 
2. 51 19E-01 
2.6607E-01 

2. 8184E-01 
2. 9854E-01 

3. 1623E-01 
3. 3497E-01 
3. 54 8 IE- 0 1 
3. 7584E-01 
3. 98 1 1E- 0 1 
4. 2170E-01 
4. 4668E-0 1 
4.7315E-01 
5. 01 19E-0 1 
5. 3088E-0 1 
5. 6234E-0 1 
5. 9566E-0 1 
6. 3096E-01 
6 . 6834E-0 1 
7. 07 94E-01 
7. 4989E-01 
7. 9433E-01 
8.4139E-01 
8. 9125E-01 
9. 4406E-01 
1.0000E 00 
1.0593E 00 
1. 1220E 00 
1 . 1885E 00 
1.2589E 00 
1.3335E 00 
1.4125E 00 
1.4962E 00 
1.5849E 00 
1.6788E 00 
1.7783E 00 
1.8836E 00 
1.9953E 00 



PHIMAX 
9.3579E-03 
1. 0500E-02 
1. 17 8 IE- 02 
1.321 8E-02 
1 . 48 3 IE- 02 
1.6641E-02 
1 . 867 IE- 02 
2. 09 50E-02 
2.35 C6E- 02 
2.63 74E-02 
2. 9592E-02 
3.3203E-02 
3. 7254E- 02 
4. 1800E-02 
4.6900E-02 
5.26 23 E-02 
5. 9044E-02 

6.62 49E-02 
7.43 32E-02 
8 . 34 02E-02 
9.3579E-02 
1.0500E-01 
1. 17 8 IE-01 
1.32 18E-01 
1.4831E-01 
1.6641E-01 
1. 8671E-01 
2. 0950E-0 1 
2. 3506E-01 

2 . 63 74 E-0 1 
2.9592E-01 
3. 32 03 E-0 1 
3.7254E-01 
4. 1800E-0 1 
4. 6900E-0 1 
5.26 23 E-0 1 
5. 9044E-01 
6.6248E-01 
7.4332E-01 
8.3402E-01 
9.3578E-01 
1.0500E 00 
1 . 178 IE 00 
1.3218E 00 
1.4831E 00 
1.6641E 00 
1.8671E 00 
2.0949E 00 
2.35C6E 00 
2.6374E 00 
2.9592E 00 
3.3203E 00 
3.7254E 00 
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WAVESLOPE SPECTRUM OF THE OCEAN 
SIGNIFICANT HEIGHT =3.3 FEET 
******** (SEA STATE 3)******** 



w 




PHIMAX 


2. 1135E 


00 


4. 1800E 


00 


2. 2387E 


00 


4. 690 OE 


00 


2.3714E 


00 


5.2623E 


00 


2.5119E 


00 


5.9044E 


00 


2.6607E 


00 


6.6248E 


00 


2. 8184E 


00 


7. 4331E 


00 


2.9854E 


00 


8. 185 1 E 


00 


3. 1623E 


00 


8. 18 5 IE 


00 


3.3496E 


00 


8. 1851 E 


00 


3. 548 1 E 


00 


8. 18 5 IE 


00 


3.7584E 


00 


8. 1851E 


00 


3. 98 1 1E 


00 


8. 18 5 1 E 


00 


4.2169E 


00 


8.1851E 


00 


4. 4668E 


00 


8. 185 IE 


00 


4.7315E 


00 


8.1851E 


00 


5.01 18E 


00 


8. 18 5 1 E 


00 


5.3088E 


00 


8. 1851E 


00 


5. 6234E 


00 


8. 18 5 1 E 


CO 


5.9566E 


00 


8. 1851E 


00 


6.3095E 


00 


8. 185 IE 


00 
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W AVESLD PE SPECTRUM OF THE OCEAN 
SIGNIFICANT HEIGHT =6.9 FEET 
********(SE£ STATE 4 ) ******** 



W 

1. 0000E-01 
1. 0593E-01 

1. 1220E-0 1 
1 . 1885E-01 
1. 2589E-01 
1 .3335E-01 
1. 4125E-01 
1 . 4962E-01 
1 . 5849 E- 0 1 
1 .6788E-01 
1. 7783E-01 
1 . 8836E-0 1 
1. 9953E-01 
2. 1135E-01 
2. 2387E-01 
2.37 14E-01 
2. 51 19E-0 1 
2. 6607E-01 
2. 8184E-0 1 

2. 9854E-01 

3. 1623E-01 
3.3497E-01 
3. 5481E-01 
3. 7584E-01 
3. 981 1E-0 1 
4.21 70E-0 1 
4. 4668E-01 
4.7315E-01 
5. 01 19E-0 1 
5. 3088E-0 1 
5. 6234E-0 1 
5. 9566E-0 1 
6. 3096E-01 
6 . 6834E-0 1 
7. 0794E-0 1 
7. 4989E-01 
7. 9433E-01 
8.41 39E-0 1 
8. 9125E-01 
9. 4406E-01 
1.0000E 00 
1.0593E 00 

1. 1220E 00 
1.1885E 00 
1.2589E 00 
1.3335E 00 
1.4125E 00 
1.4962E 00 
1.5849E 00 
1.6788E 00 
1.7783E 00 
1.8836E 00 
1.9953E 00 



PHIHAX 
1 . 9566E-02 

2. 19 54 E-02 

2. 46 33E-02 
2. 7638E-02 
3. 10 1 1E-02 

3 . 47 95 E-02 
3. 9040 E-02 
4.38 04 E-02 
4. 91 49E-02 
5.5146E-02 
6. 1875E- 02 
6. 9424E-02 
7.78 S5E- 02 
8.74 00 E-02 
9. 8065E-02 
1. 1003E-01 
1. 2346E-01 
1 .3852E-01 
1 . 5542E- 0 1 
1.7439E-01 
1 . 9566E-0 1 
2. 1954E-01 
2.46 33E-0 1 
2.7638E-01 

3. 101 IE-01 
3. 4795E-01 
3. 9040E- 01 
4.3804 E-0 1 

4. 9149E-01 
5. 51 46E-01 
6. 1874E-C1 
6.94 24 E-0 1 
7.78 S5E- 0 1 
8 . 74 00 E-0 1 
9. 8064E-0 1 
1.1003E 00 
1.2346E 00 
1.3852E 00 
1.5542E 00 
1.7439E 00 
1.9566E 00 
2. 19 54 E 00 
2.4633E CO 
2.7638E 00 
3.1011E 00 
3.4794E 00 
3.9040E 00 
4.3803E 00 
4.9148E 00 
5. 51 45E 00 
6. 1874E 00 
6.9424E 00 
7.78S5E 00 



\ 
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WAVESLOPE SPECTRUM OF THE OCEAN 
SIGNIFICANT HEIGHT =6.9 FEET 
******** (SEh STATE 4) ******** 



w 




PHIMAX 


2. 1135E 


00 


8. 1851E 


00 


2. 2387E 


00 


8.1851E 


00 


2.3714E 


00 


8. 1 8 5 1 E 


00 


2. 51 19E 


00 


8. 185 IE 


00 


2.6607E 


00 


8. 1851E 


00 


2. 8184E 


00 


8. 1851E 


00 


2. 9854E 


00 


8. 1851 E 


00 


3. 1623E 


00 


8. 18 5 1 E 


00 


3.3496E 


00 


8. 185 1 E 


00 


3. 548 IE 


00 


8. 18 5 1 E 


CO 


3. 7584E 


00 


8. 1851E 


00 


3. 98 1 1E 


00 


8. 18 5 1 E 


00 


4.2169E 


00 


8. 1851 E 


00 


4. 4668E 


00 


8. 1851E 


00 


4.7315E 


00 


8. 185 1 E 


00 


5.01 18E 


00 


8. 185 1 E 


00 


5.3088E 


00 


8. 185 IE 


00 


5.6234E 


00 


8. 1851E 


00 


5.9566E 


00 


8. 185 IE 


00 


6.3095E 


00 


8. 185 IE 


CO 
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WAVESLOPE SPECTRUM OF THE OCEAN 
SIGNIFICANT HEIGHT = 10. FEET 
******** (sea STATE 5) ******** 



W 

1. 00 00E-01 
1 . 0593E-0 1 

1. 1220E-01 
1 . 1 885E-0 1 
1. 258 9E-0 1 
1 .3335E-01 
1. 4125E-01 
1.4962E-01 

1. 5849E- 0 1 
1 .6788E-01 
1. 77 83E-01 
1 . 8836E-0 1 
1.9953E-01 

2. 1135E-01 
2. 2387E-01 
2. 3714E-01 
2. 5119E-01 
2. 6607E-01 
2. 8184E-01 
2. 9854E-01 

3. 1623E-01 
3.3497E-01 
3. 548 IE- 0 1 
3 . 7584E- 0 1 
3. 9811E-01 
4.217 0E- 0 1 
4. 4668E-0 1 
4.7315E-01 
5. 01 19E-0 1 

i 5 . 3088E- 0 1 

5. 6234E-0 1 
5. 9566E-0 1 

6. 30 96E-01 
6. 6834E-01 
7. 0794E-0 1 
7. 4989E-01 
7. 9433E-0 1 
8. 4139E-01 
8.91 25E-0 1 
9. 4406E-01 
1.0000E 00 
1.0593E 00 

1. 1220E 00 
1.1885E 00 
1.2589E 00 
1.3335E 00 
1.4125E 00 
1.4962E 00 
1.5849E 00 
1.6788E 00 
1.7783E 00 
1.8836E 00 
1.9953E 00 



2. 8357E-02 

3. 18 17E-02 
3. 5700E- 02 
4. 0056E-02 
4. 4943E-02 
5. 0427E-02 
5.6580E-02 
6. 3484E-02 
7. 1230E-02 
7.992 IE-02 
8. 9673E-02 
1. 0062E-01 
1. 1289E-01 
1.2667E-01 
1. 4212E-01 
1.5946E-01 
1.7892E-01 
2. 0075E-01 

2. 25 25E-01 
2. 52 73E-01 

2. 8357E-01 

3. 1817E-01 
3. 5699E-01 
4. 0056E-01 
4. 4943E- 0 1 
5. 0427E-0 1 
5. 6580E-01 
6.3484E-01 
7. 1230E-01 
7. 9921E-01 
8. 9673E-01 
1.0061E 00 
1.1289E 00 
1.2667E 00 
1.4212E 00 
1.5946E 00 
1.7892E 00 
2.0075E 00 
2.2525E 00 
2.5273E 00 
2.8357E 00 
3.1817E 00 
3.5699E 00 
4.0055E 00 
4.4943E 00 
5.0427E 00 
5.6580E 00 
6.3483E 00 
7.1230E 00 
7. 9921 E 00 
8.1851E 00 
8. 1851 E 00 
8.1851E 00 
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WAVESLOPE SPECTRUM OF THE OCEAN 
SIGNIFICANT HEIGHT = 10. FEET 
****###*(SEA STATE 5) *#**#**# 



W 

2.1135E 00 
2.2387E 00 
2.3714E 00 
2.5119E 00 
2.6607E 00 
2.8184E 00 
2.9854E 00 
3. 1623E 00 
3.3496E 00 
3.5481E 00 
3.7584E 00 
3.9811E 00 
4.2169E 00 
4.4668E 00 
4.7315E 00 
5. 01 18E 00 
5.3088E 00 
5.6234E 00 
5.9566E 00 
6.3095E 00 



PHIMAX 
8.1851E 00 
8.1851E 00 
8.1851E 00 

8. 185 IE 00 
8. 185 1 E 00 
8.1851E CO 
8.1851E 00 
8. 18 5 1 E 00 
8 . 185 1 E 00 
8.1851E 00 
8. 18 5 1 E 00 

8. 185 IE 00 
8. 18 5 1 E 00 
8. 185 IE 00 
8. 1851 E 00 
8. 185 IE 00 
8. 1851E 00 
8.1851E CO 
8.1851E 00 

8 . 185 IE 00 



78 



WAVESL3PE SPECTRUM OF THE OCEAN 
SIGNIFICANT HEIGHT = 15. FEET 
********( SEA STATE 6) ******** 



W 

1. OOOOE-O 1 
1 . 0593E-01 
1. 1220E-01 
1. 1885E-01 
1. 2589E-01 
1 . 3335E-0 1 
1.4125E-01 
1.4962E-01 
1. 5849E-01 
1 .6788E-01 
1. 7783E-01 
1 . 8836E-0 1 
1. 9953E-01 
2. 1135E-01 
2. 2387E-01 
2. 3714E-01 

2. 51 19E-01 
2 . 6607E-0 1 
2. 8184E-01 
2. 9854E-01 

3. 1623E-01 
3.3497E-01 
3. 54 8 IE- 01 
3. 7584E-01 
3. 9811E-01 
4.21 70E-0 1 
4.46 68E- 0 1 
4.731 5E-0 1 
5. 0119 E- 01 
5 . 3088E-0 1 
5. 6234E- 0 1 
5. 9566E-01 

6. 3096E-01 
6 . 6834E-0 1 
7. 0794E-01 
7 . 4989E-0 1 

7. 9433E-0 1 
8. 4139E-01 
8. 9125E-0 1 
9. 4406E-01 



1 , 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 , 



OOOOE 
0593E 
1220E 
1885E 
2589E 
3335E 
41 25E 
4962E 
5849E 
6788E 
77 83E 
8836E 
9953E 



00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 



FHIMAX 
4. 2536E-02 
4. 7726E-02 
5.3549E-02 
6 . 0083E-02 
6.74 15E-02 
7.5641E-02 
8.4870E-02 
9.5226E-02 
1.0685E-01 
1 . 1988E-01 
1 . 34 5 IE- 0 1 
1 . 50 92E-0 1 
1. 6934E-01 
1.90 OOE-O 1 
2. 1318E-01 
2.3920E-0 1 
2. 6838E-01 
3. 01 13E-0 1 
3. 3787E-01 
3.79 10E-01 
4. 2536E-0 1 
4.7726E-01 
5. 3549E-01 
6. 0083E-0 1 
6.74 14E-01 
7. 5640E-01 
8. 4870E-01 
9.5225E-01 
1.0684E CO 
1.1988E 00 
1.3451E 00 
1.5092E 00 
1.6934E 00 
1.9000E 00 
2. 13 18E 00 
2.3920E 00 
2.6838E 00 
3.0113E 00 
3.3787E CO 
3.7910E 00 
4.2536E 00 
4.7726E 00 
5.3549E 00 
6.0083E 00 
6.7414E 00 
7.5640E 00 
8. 185 IE 00 
8. 1851 E 00 
8.1851E 00 
8. 1851E 00 
8. 18 5 IE 00 
8. 185 1 E 00 
8. 185 IE 00 
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WAVESLOPE SPECTRUM OF TEE OCEAN 
SIGNIFICANT HEIGHT = 15. FEET 
*#*#***« (SEA STATE 6 ) ******** 



W PHIMAX 



2.1135E 00 
2. 2387E 00 
2.3714E 00 
2.5119E 00 
2.6607E 00 
2.8184E 00 
2.9854E 00 
3. 1623E 00 
3.3496E 00 
3.5481E 00 
3.7584E 00 
3. 98 1 1E 00 
4.2169E 00 
4.4668E 00 
4.7315E 00 
5.0118E 00 
5.3088E 00 
5.6234E 00 
5.9566E 00 
6.3095E 00 



8. 


185 


IE 


00 


8. 


185 


IE 


CO 


8. 


185 


IE 


00 


8. 


185 


IE 


00 


8. 


185 


IE 


00 


8. 


185 


IE 


00 


8. 


1851E 


00 


8. 


185 


IE 


00 


8. 


185 


IE 


00 


8. 


185 


IE 


CO 


8. 


185 


IE 


00 


8. 


185 


IE 


00 


8. 


185 


IE 


00 


8. 


1851E 


CO 


8. 


185 


IE 


00 


8. 


185 


IE 


00 


8. 


185 


IE 


00 


8. 


185 


IE 


00 


8. 


185 


IE 


00 


8. 


185 


IE 


CO 
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